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*vj \ Abstract. We determine the evolution from z = 3 — s- of the galaxy 

f^ ■ and mass correlation functions and bias factor in a 50/i~^Mpc ACDM hy- 

drodynamic simulation with 10/i~^kpc resolution. The mass correlation 
^^ • function grows with time, but the galaxy correlation function shows little 

evolution and is well described by a power law. At early times, galaxies 
,-S^ ■ are biased traces of mass, with bias being higher on smaller scales. By 

^i z = 0, galaxies trace the mass, and the bias shows little scale dependence. 
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'^ I The correlation function of galaxies, ^g, is a standard observational mea- 

Ci ■ sure of structure in the Universe. Understanding the relationship between this 

observable and the correlation of the underlying mass distribution, ^m, is cru- 
cial for relating observations of galaxy clustering to predictions of cosmological 
rS I models. The relationship between ^g and ^m is usually characterized through 

j^ ■ the bias parameter b = ag/am, where a is the rms fluctuation on some specified 

scale. Here we investigate the evolution of ^g, ^m and 6 in a ACDM simulation. 
We simulate a random 50/i~^Mpc cube in a ACDM universe, with Q.m = 0.4, 
Qa = 0.6, Hq = 65, n = 0.95, and Qfe = 0.02/i-2. We use Parallel TreeSPH 
to advance 144^ gas and 144^ dark matter particles from z = 49 to z = 0. 
Our spatial resolution is 10/i~^ kpc, and our mass resolution per particle is 
m-SPH = 8.5 X 10^ Mq and rudark = 6.3 x W^Mq. Using a 60-particle criterion 
for our simulated galaxy completeness limit [1] implies that we are resolving most 
galaxies with Mi^aryonic ^ 5 X 10^^ Mq. We include star formation and thermal 
feedback. Galaxies are identified using Spline Kernel Interpolative DENMAX. 
At 2 = 3, 2, 1, and 0, we identify 929, 1781, 4138, and 5264 galaxies, respectively. 
Note that the majority of galaxies form between z = 2 and z = 1 in this model. 
We consider the 500 most massive galaxies at each redshift, representing 
a constant comoving galaxy density of n = 0.004/i'^Mpc~^, comparable to that 
of Lyman break galaxies and present-day L* galaxies. The galaxy correlation 
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function ^g at each redshift is shown as open circles in Figure 1. A power-law 
fit to these points is shown as the solid line, and the best-fit values for the 
correlation length tq and slope 7 are indicated in the upper right. The mass 
correlation function £^rn is shown as the dashed line. 
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Figure 1 (left 4 panels): Correlation function of galaxies and mass at various redshifts. 

Figure 2 (right panel): Galaxy bias on various scales as a function of redshift. 

^g evolves little from z = 3 to z = 0, in agreement with previous studies [2]. 
tq is roughly constant, with a slight increase to lower redshifts. 7 also evolves 
only slightly, with ^g becoming steeper at low redshift. Conversely, ^m increases 
substantially with time. The z = correlation function of this model is in good 
agreement with observations. 

The difference in evolution of S^g and S,m is reflected in the evolution of 
b{R). We consider 4 scales, i? = 2,4, 8 and 12/i~^Mpc. The evolution of the 
bias parameter on these scales is shown in Figure 2. At early times, galaxies 
are highly biased tracers of the underlying mass distribution, in agreement with 
previous studies(see e.g. [3]). This implies that Lyman break galaxies do not 
trace the underlying mass distribution at high z. For R ^ 8/i~^Mpc, h increases 
to smaller scales at high z. The bias factor declines with time on all scales, and 
shows little scale dependence for z <^ 1. By z = 0, 6 ~ 1.2 on all scales from 
R = 2 —>■ 12/i~^Mpc, reflecting cg = 0.8 chosen for our simulation. 
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